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ABSTRACT 

Context. Determining the mass transfer in a close binary system is of prime importance for understanding its evolution. SS Leporis, 
a symbiotic star showing the Algol paradox and presenting clear evidence of ongoing mass transfer, in which the donor has been 
thought to fill its Roche lobe, is a target particularly suited to this kind of study. 

Aims. Since previous spectroscopic and interferometric observations have not been able to fully constrain the system morphology 
and characteristics, we go one step further to determine its orbital parameters, for which we need new interferometric observations 
directly probing the inner parts of the system with a much higher number of spatial frequencies. 

Methods. We use data obtained at eight different epochs with the VLTI instruments AMBER and PIONIER in the H- and K-bands. 
We performed aperture synthesis imaging to obtain the first model-independent view of this system. We then modelled it as a binary 
(whose giant is spatially resolved) that is surrounded by a circumbinary disc. 

Results. Combining these interferometric measurements with previous radial velocities, we fully constrain the orbit of the system. 
We then determine the mass of each star and significantly revise the mass ratio. The M giant also appears to be almost twice smaller 
than previously thought. Additionally, the low spectral resolution of the data allows the flux of both stars and of the dusty disc to be 
determined along the H and K bands, and thereby extracting their temperatures. 

Conclusions. We find that the M giant actually does not stricto sensus fill its Roche lobe. The mass transfer is more likely to occur 
through the accretion of an important part of the giant wind. We finally rise the possibility for an enhanced mass loss from the giant, 
and we show that an accretion disc should have formed around the A star. 



Key words. Techniques: interferometric - Binaries: symbiotic 
mental parameters - Accretion, accretion disc 
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1. Introduction 

Symbiotic stars are interacting binaries composed of a hot star 
accreting material from a more evolved red giant companion. 
They are excellent laboratories for studying a wide specttum of 
poorly understood physical processes, like the late stage of stel- 
lar evolution, the mass loss of red giants, and the mass trans- 
fer and accretion in binary systems (Mikolaiewska 2007). Their 
study has important implications for a wide range of objects, 
like Type la supernovae, barium stars, the shaping of plane- 
tary nebulae, and compact binarie s like cataclysmic variables 
dPodsiadlowski & Mohamedl2007l) . 

SS Leporis (17 Lep; HD 41511; HR 2148) is a prime exam- 
ple of such a long-period interacting system, even though it does 
not belong to the most common symbiotic systems, because the 
hot star is not the usual compact white dwarf. As such, SS Lep 
is a symbiotic system in the first phase of mass transfer, while 
most symbiotic stars are in their second episode of mass trans- 
fer, following the first one that produced the white dwarf. 

SS Lep has been known for many decades to present 
symbiotic features, and its optical spectrum shows at least 



* Based on observations made with the VLTI European Southern 
Observatory telescopes obtained from the ESO/ST-ECF Science 
Archive Facility. 



three components Jsiruve & Roachl Il939t lMolaroetaflll983b 
Weltv & Wade 19951). The spectral lines of an A star are largely 
obliterated by shell features that dominate at shorter wave- 
lengths, while an M star spectrum becomes inc reasingly obvious 
at longer wavelengths. Welt v & Wadd (1 19951) estimated an M4 
III spectral type for the cool companion, while even earlier types 
have been estimated by previous authors. The shell is absorbing 
light primarily from the A star, indicating some mass loss from 
the hotter star. The system, however, presents the so-called Algol 
paradox, as the most evolved star is also the least massive, which 
implies that the cool star must have lost a large quantity of mat- 
ter and that the hot companion ha s accreted part or most of it. 
More over, the regular "outb ursts" dStruvdl 19301: Weltv & Wade 
119951) and the UV activity dPolidan & Shordl 1993b of the A star 
shell are clear testimony t on going mass-transfer e pisodes. From 
interferometric observations. lVerhoelst et afl(l2007l) inferred that 
the mass transfer occurs because the M giant fills its Roche lobe. 

The binary system is additionally surro unded by a larg e 
circumbinary dust disc and/or envelope dJura et alJ [2001 ). 
Interferometric observations confirmed t his fact by revealing it s 
presence in the inner part of the system (IVerhoelst et al.ll201 ll) . 
further noticing that the structure must be in a disc-like geome- 
try to be comp a tible w ith the low extinction towards the central 
star. lJura et all (12001 1) suggest that the circumbinary disc con- 
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Table 1. Previously estimated parameters of SS Lep 



System 




A star 


M star 


d [pc] 


279 ± 24 " 


SpT 


A1V 


M6III 1 


P[d] 


260.3 ± 1.8* 


T cB [K] 


9000 c4 


3500 c - d 


e 


0.024 + 0.005 b 


9 [mas] 


0.53 ± 0.02 c 


3.11 +0.32 c 


i 


30° + 10° c 




-15 c 


110 ± 30 c 


f(M) 


0.261 ± 0.005 >' 


F [%] 


1 1 ± 7 c 


60 ± 50 c 


l/q 


4 ± 1 *- c 


M[M Q ] 


2-3 a - b 


0.4-1 "•*> 



Notes, d is the distance, P the orbital period, e the eccentricity, ! the 
inclination, f(M) the mass function, q = M M /M A the mass ratio. For 
the stars, SpT is the spectral type, r e g the temperature, the apparent 
diameter, R the linear radi us, F the flux cont ribu tion at 2.2 fim, and M 
the mass. References: (a ) IVan Leeuwerill2007l: (b ) IWeltv & Wadell 19951 : 
(c) IVerhoelst et alj2007l : (dl lBlondel et al J 19931 . 

tains large grains that are formed by coagulation and, based on 
the large and rather unique 12 jjm excess of SS Lep despite its 
rather low luminosity, that the disc may be losing mass by a wind 
at a rate of 8 x 1(T 9 M yr~. 

The orbital characteristics and circumbinary disc of SS Lep 
very closely resemble those of the pos t-AGB binaries wit h sta- 
ble, Keplerian circumbinary dust discs (Ivan Winckell2 003). and, 
as such, SS Lep may be considered as a system linking binary M 
giants and post-AGB systems - the M-giant should indeed very 
soon evolve into a post-AGB star. In those post-AGB binaries, 
the spectra are rich in crystalline features while the spectrum 
of SS Lep appears entirely amorphous, which, if a link is in- 
deed in order, would imply further dust processing in the disc. 
Whether the disc can survive long enough if it is indeed losing 
mass through a win d is, however, still an open question. 

ICowlevi (1 19671) argued that the system consists of a B9 V 
primary and an Ml III secondary in a 260-d eccentric orbit 
(e = 0.132). She developed a scenario in which the secondary 
fills its Roche lobe near periastron and mass transfer proceeds 
for a short time ther eafter. From spectra covering 3.5 orbits, 
IWeltv & Wadel (1 19951) proved this scenario unlikely, as their re- 
vised orbit provided a similar orbital period but a significantly 
reduced eccentricity e = 0.024. They estimated a mass ratio of 
l/q = 3.50 + 0.57, where the error was very likely severely un- 
derestimated given the poor fit of the single Mg II line they used 
to measure th e radial velocity of th e A star. 

Recently, IVan Leeuwenl (120071) has reevaluated the parallax 
of SS Lep from Hipparcos data, obtaining n - 3.59 + 0.31 mas, 
that is, a distance of 279 ± 24 pc, so smaller than the previously 
and generally used value of 330^ pc. The most recent param- 
eters of SS Lep as collected from the literature until the present 
work, is presented in Table [T] 

We report here interferometric observations that have al- 
lowed model-independent image synthesis (Section [3]) and a 
more precise modelling of the system (Section|4]i. Section[5]fo- 
cuses on the binary by determining mainly the orbital parame- 
ters and the M star diameter. We finally discuss the mass transfer 
process in Section|6] 

2. Observations 

Data w ere collected at the Very Large Telescope Interferometer 
VLTLlHaguenauer et al]|2010l) with the spectrograph AMBER 
Petrov et alJ 12007]) and the four- telescope visitor instrument 
PIONIER (iLe Bouquin et al.ll201 ll) . All observations made use 
of the 1.8-m Auxiliary Telescopes. Table |2]presents the observa- 
tion log. Three data sets are the combination of two observations 



Table 2. VLTI Observation log of SS Lep for the AMBER and 
PIONIER observations. 



Run 


Date 


Baselines 


Range 


Al 


11-11-2008 


E0-G0-H0 


15 - 


130m 




13-11-2008 


A0-G1-K0 






A2 


26-12-2008 


A0-G1-K0 


90- 


130 m 


A3 


21-02-2009 


D0-G1-H0 


15 - 


75 m 




28-02-2009 


E0-G0-H0 






A4 


07-04-2009 


D0-G1-H0 


65 - 


70 m 


PI 


28-10-2010 


D0-G1-H0-I1 


15 - 


80 m 




30-10-2010 


D0-E0-H0-I1 






P2 


29-11-2010 


E0-G0-H0-I1 


15 - 


70 m 


P3 


07-12-2010 


D0-G1-H0-I1 


45 - 


80 m 


P4 


22-12-2010 


A0-G1-K0-I1 


45 - 


130m 



separated in time by a few days (seven to the maximum) to in- 
crease the (M,v)-plane coverage. We made sure that this brings 
valuable information without biasing the results, more specif- 
ically the estimation of the binary orientation. Typical (u,v)- 
planes for AMBER and PIONIER observations can be seen in 
Fig.ffl 

AMBER data We used archive AMBER data obtained during 
four different nights in a period of 200 days (more than half an 
orbital period). They cover simultaneoulsy the J-, H-, and K- 
bands with a spectral resolution R ~ 35. Even though J-band 
fringes have been properly recorded in several observations, we 
decided to discard them from the analysis of this paper since 
the data quality is significantly worse than for longer wave- 
lengths. Raw visibility and closure phase values were computed 
using the latest pub lic version of the amdlib package (version 
3; iMalbet et al.ll2.Q10l) and the yorick interface provided by the 
Jean-Marie Mariotti Center. 

PIONIER data They were obtained in the H-band between the 
end of October 2010 and December 2010 during the commis- 
sioning runs of the instrument. We used the prism that provides 
a spectral resolution R ~ 40, that is, six spectral channels across 
the H-band. Because these observations have been made during 
the commissioning, the (M,v)-plane coverage is still relatively 
poor for a whole night of observations but more complete 
than AMBER' s. Data were reduced with the pndrs package 
presented bv lLe Bou quin et al] d2011l) . 

The low spectral resolution multiplies the number of spatial 
frequencies, and brings a wealth of information (Sect. [4]). This is 




-50 50 -50 50 



U[10 6 rad-'] U[10 6 racr'] 

Fig. 1. Typical («,v)-plane coverage for AMBER (left, 
November 2008) and PIONIER (right, October 2010) ob- 
servations. 
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28-10-2010 07-12-2010 



22-12-2010 



- 




A star 



M giant 




-5 



E (mas) 



E (mas) 



E (mas) 



Fig. 2. Model-independent image reconstruction of SS Lep obtained during the PIONIER runs PI, P2, and P4. The resolved M giant 
and the A star are clearly identified. The images are centered on the center of mass (central cross) as determined from Section IBT2l 
The distortion of the giant in the image is most certainly due to an asymmetric PSF rather than to a definite tidal effect. Three faint 
artefacts are visible on the periphery of the image. 



especially true with AMBER for which the (M,vJ-coverages were 
relatively poor. 



3. Image synthesis 

With its four telescopes, PIONIER provides six visibilities and 
four closure phases simultaneously, which allows a reliable im- 
age reconstruc tion for the fou r observations. We used the MIRA 
software from lThiebauj d2008l) . MIRA proceeds by direct min- 
imisation of a penalised likelihood. This penalty is the sum of 
two terms: a likelihood term that enforces agreement of the 
model with the data (visibilities and closure phases), plus a reg- 
ularisation term to account for priors. The priors are required to 
lever the possible degeneracies due to the sparseness of the spa- 
tial frequency sampling. We use here the "total variation" regu- 
lar isation associated w ith positivity constraint as recommended 
bv lRenard et al.l(l201 ll) . The pixel scale is 0.25 mas and the field- 
of-view is 200x200 pixels. The starting point is a Dirac function 
in (0,0). We set the hyper-parameter to a low value of 100, so 
that the weight of the regularisation term is kept small with re- 
spect to the fit to the data. It brings some super-resolution, at the 
cost of an improved noise level in the image. We combined all 
the spectral channels to improve the (u, v)-plane coverage. That 
the image is indeed "grey" over the H-band is demonstrated in 
the next section. The reconstructed images for runs PI, P2, and 
P4 are presented in Fig. [2] Each image shows the binary nature of 
SS Lep, the separation being slightly smaller than 5 mas. From 
one observation to the next, we can observe the rotation of the 
system. 

The A st ar and its shell have an expected spatial extension 
of 0.5 mas dVerhoelst et al.ll2007h so that we do not expect to re- 
solve them with our VLTI baselines. Therefore, the size of the 
spot corresponding to the A star more or less defines the point 
spread function (PSF) of the image, about 1 mas large. Because 
the M giant is the most luminous component of the system in 
the H-band, we identify it in the image as the darkest spot. With 
respect to the A star, we clearly see that it is spatially resolved 
and measures approximately 2 mas in diameter. We expect the 



distortion observed in the image to come from an asymmetric 
filling of the (u, v)-plane (implying a non-circular PSF on the 
reconstructed image) rather than to a real tidal distortion. As a 
matter of fact, the tidal distortion would be around 5 ~ 7%, 
i.e. less than seen in the image. Additionally, its orientation in 
the image corresponds well with the asymmetry observed in the 
corresponding (u, v)-planes. It was actually not possible to im- 
age the circumbinary disc because of the lack of data with short 
baselines. 



4. Modelling 

Our observations clearly show that SS Lep is a spatially resolved 
binary whose M giant is actually resolved for all observations. 
We built a geometrical model to determine the characteristics of 
the individual components. The M giant and the A star are mod- 
elled as uniform discs, and the circumbinary material is mod- 
elled as a Gaussian envelope. We tried to detect a possible tidal 
distortion of the giant or matter escaping from its atmosphere by 
modelling it with an elongated uniform disc. Results were not 
persuasive and, similar to the image reconstruction, we cannot 
conclude anything about this because we lack the longest base- 
lines able to measure distortion of a few percent. The spatial 
resolution of 1 mas was also not sufficient to resolve the putative 
shell or an accretion disc around the A star, which agrees with 
the 0.5 m as size estimated fro m the spectral energy distribution 
(SED) in lVerhoelst etal.1 d2007l) . We therefore fixed its diameter 
to 0.5 mas. 

The model we used to fit the interferometric data (visibilities 
and closure phases) therefore comprises six degrees of freedom: 
the relative flux contribution of two components of the system; 
the binary separation and its orientation, the size of the M giant, 
the size of the circumbinary envelope. Our data sets are perfectly 
suited to spectral analysis. To properly fit the data it appeared 
necessary to consider the fluxes to be wavelength-dependent. 

The data and the results of our fits are presented in the 
Appendix. Starting from the PIONIER images, we are able to 
measure the binary separation and orientation for each dataset 
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Table 3. Orbital paramete rs of SS Lep obtained by combining 
previous radial velocities (IWeltv& Waddll995h with our 8 as- 
trometric measurements. 



-4 -2 2 4 

Right Ascension (milli-arcsecond) 



Fig. 3. SS Lep best orbit ( dashed line) obtained by combining 
previous radial velocities (IWeltv& Wadd 119951) with our as- 
trometric measurements. The central dot indicates the A star. 
AMBER and PIONIER points are respectively presented by the 
red and blue crosses representing the 3-cr error bars. The cor- 
responding points on the best orbit are indicated by the short 
segments originating in each point. 



independently. The relative flux of the three components could 
be recovered between 1.6 and 2.5 fim, with a dispersion of 3% 
between the different epochs. The relative flux ratios are al- 
most constant over the H-band, which validates the "grey" ap- 
proach used in Sect. [3] We tried to measure a chromatic diam- 
eter for the giant but results were not consistent between the 
different epochs. Finally, despite the relatively long period be- 
tween AMBER and PIONIER observations (almost two orbital 
periods), we note rather good consistency of results within error 
bars, indicating that the system is relatively stable. 

We were not able to extract much information about the 
large circumbinary disc because of the lack of very short base- 
lines. The disc is almost totally resolved with the smallest pro- 
jected baselines of 15 m. We were only able to roughly de- 
termine its full width at half maximum (FWHM) for only 
two observations out of eight (Al and A4), and its relative 
flux for 6 of them (Al, A3, PI to P4). We measured an 
FW HM of 12.2 + 0.2 mas in agreement with the estimation of 
IVerhoelst et al.l d2007l) i n the near-IR. The di sc has been ob- 
served in the mid-IR by IVerhoelst et al.1 (1201 ll) . who measured 
a Gaussian FWHM of 26 mas. Obervations with the 10-metre 
baselines of the Keck Segment- Tiltin g Experiment at 10.7 /im 
did not resolve it ( Monni er et al.l20 09) and indicate that it should 
not be larger than 200 mas in the mid-IR. 



5. Characteristics of the individual components 

5.1. The orbit of SS Lep 

To compute the most reliable orbit possible, we combined our 
eight astrometric positions of the binary with the radial veloc- 
ities of the M star obtained by IWeltv & Wadel (fl995fl We 
deduced the orbital parameters from a global x 2 minimisation 
of these data. The best-fit orbit is shown in Fig. [3] and the or- 
bital parameters are listed in Table [3] Uncertainties on the or- 
bital elements are estimated via Monte-Carlo simulations. The 



Semi major axis a 
Linear semi major axis a 
Inclination i 
Eccentricity e 

Longitude of the ascending node CI 
Argument of periastron a> 



4.492 ±0.014 mas 
1.26 ± 0.06 AU 
143.7 ± 0.5° 
0.005 ± 0.003 
162.2 ± 0.7° 
118 ±30° 



inclination angle of 143.7° and the non -significant ecc entricity^ 
confirm the previous measurements of Welty & Wade] (I 1995b of 
a quasi-circularised orbit observed with an inclination angle be- 
tween 28° and 38°. The circularised orbit is not a surprise for 
an evolved sy mbiotic system, with an M giant in a short orbit. 
For instance. lFekel et al. I d2007l) find that 17 of the 21 (i.e. 81%) 
red symbiotic systems with periods P < 800 days have circular 
orbits. 

This result defi nitely inv alidates the periastron-passage mass 
transfer scenario o f lCowlevl(fl967l) . which required a significant 
eccentricity (e = 0.134) to explain the regular "outbursts" of 
the system. Finally, given the almost null eccentricity, co is also 
poorly constrained. 



5.2. The masses 

Combining our valu e for the inclination w ith the binary mass 
function obtained bv lWeltv& Wad3(fT99l . we can estimate the 
individual mass of the stars, hence the mass ratio. The main 
source of uncertainty in this estimation resides in the distance, 
as determined by Hipparcos. Using the distance and the angu- 
lar separation of the two stars, we obtain a = 1.26 + 0.06 AU, 
and thus, through Kepler's third law, the total mass of the sys- 
tem is estimated as 4.01 + 0.60 M . We then derive Ma = 
2.71 + 0.27 M Q , M m = 1.30 + 0.33 M Q , and 1/q = M A /M M = 
2.17 + 0.35. The mass ratio is thus much greater than previously 
thought. While the A star still has a mass in the range 2-3 M Q , 
the mass of the M giant is now much higher than estimated ear- 
lier, and this implies that less matter was transferred into the 
system than previously guessed. We come back to this later. 

5.3. The M star 

Averaging over all the epochs, we measure an apparent diameter 
for the M star Om.ud - 2.208 + 0.012 mas, where the error is 
computed from the dispersion of the eight estimations. The error 
bars do not include systematic effects - e.g. due to any tidal dis- 
tortions -, and could be a few percent. It was also not possible to 
identify a dependence in the giant size as a fu nction of the wa ve- 
length. The previous VINCI observations of Verhoelst et aT]led 
to a higher value of 2. 94+0. 3 mas, most likely because theirs was 
the result of a one-year survey of the source, without any phase 
information in the interferometric data. This involved modelling 
the system as a symmetric object, so that it was impossible to dis- 
entangle the signatures of the rotating binary from the resolved 
giant one. 

The conversion factor from the uniform disc to a limb- 
darkened one differs by a few percent depending on the au- 



We discarded the A star radial velocitites that were not convincing. 



2 The astrometric points alone lead to an eccentricity of 0.004±0.008, 
compatible with a circular orbit. 
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1 ' 1 1 1 Table 4. Stellar parameters extracted in this study. 




" 5 ; w+mt-H + 




X [microns] 

Fig. 4. Flux of the M giant (red), the A star (blue), and the en- 
velope (magenta). The grey curve is the M star MARCS spec- 
trum. In black is the sum of the three components adjusted to the 
2MASS magnitudes in the H- and K-bands. The dots are the data 
plus the error bars, and the solid lines are the models for each of 
the components. 



thorfl We adopt a value of 1.04, which leads to a limb-darkened 
diameter equal to 9m,ld = 2.296 + 0.013 mas. These results 
ag ree with the l imb-darkened diameter estimated from the SED 
of IVerhoelst et alj (6m,ld = 2.66 + 0.33 mas). Taking the un- 
certainty on the distance into account, the M giant radius is 
Rm = 66.7 + 3.3 R Q , which is 40% smaller than previously ob- 
tained. This leads to a surface gravity logg ~ 0.9. 

iDumm & Schildl (1 19981) provide measurements of radius for 
(non-Mira) M star, showing that with a radius around ~ 65 R , 
the M star spectral type shoul d be M1/M2. Given the orbital pe- 
riod. iMtirset & Schmidl d 1999) agree by deriving a spectral type 
between M0 and Ml. From the table of these authors, stars with 
similar luminosities to the one estimated for the M giant of SS 
Lep all have a radius between 42 and 67 Rq. 

Based on the dependence of the nuclear time scale on stel- 
lar mass, the M star must have an initial mass at least 20% 
higher than its companion to have evolved on the AGB, while 
its companion is still in the phase of central hydrogen burn- 
ing. Because the initial mass of the system was probably greater 
than its current value, this implies an initial mass of the giant 
Mm,o > 2.2 M 0. Finally, stars with such masses (unless much 
larger) never go through a stage with large radii when on the 
RGB - with the maximum radius reached of the order of 30 R 
- which indicates the M star is more likely on the early-AGB 
phase. 





M star 


A star 


Mass [Mo] 


1.30 ± 


0.33 


2.71 ±0.27 


Apparent diameter [mas] 


2.208 ± 


0.012 


0.6 ± 0.05 


Linear radius [R Q ] 


66.7 ± 3.3 


~ 18 


Temperature [K] 


3500 d 


i 200 


-9000 



5.4. Temperature of individual components 

The modelling presented in Sect. [4] allows us to estimate the 
relative flux of the three components between 1.6 and 2.3 yum 
(see Appendix). These measurements can be used to constrain 
the temperature and the size of the individual components. 
We impose that the M giant diameter tiso 2.2 mas as esti- 
mated previously. We model its SED with a MARCS model 
whose dependencies ar e the temperature and the metallicity 
dGustafsson et al.1 120 08). The A star is in the Ray leigh- Jeans 
regime: it is impossible to fit simultaneously its temperature 
and size. We leave its diam eter free and model its SEP by a 
KURUCZ model at 9000 K dCasteUi & Kuruczll200l . The disc 
is modelled by a Gaussian whose FWHM and bl ackbody tem- 
peratu re are left free. We use the absorption law oflCardelli et aTl 
(119891) with R v = 3.1 and A v = 0.7 mag (IVerhoelst et al.l l2007; 
Malfait, Bogae rt. & Waelkenslll998l) . Additionally, we force the 
total SED (M giant + A st ar + disc) to be co mpatible with 
the 2MASS measurements (ISkrutskie et alj|2006l) . The best fit 
is presented in Fig. [4] The lack of an absolute spectrum makes 
difficult a more realistic modelling. 

The M star temperature is found to be around 3500 + 200 K. 
We also confirm that the A star is apparently larger than expected 
from its spectral class (9a - 0.6 + 0.05 mas, or a linear radius of 
about 18 R Q ). For the disc, we found a blackbody temperature 
of 1700 ± 100 K and an FWHM of 8.0 + 0.5 mas. Interestingly, 
this is incompatible with the 12 mas derived from the fit of vis- 
ibility curves in Sect. [4] We see this inconsistency as a hint that 
a Gaussian geometry is probably too simple to model the cir- 
cumbinary environment. 

6. The mass transfer process 

As explained above, SS Lep shows evidence of mass transfer be- 
tween the M giant and the A star, and this mass transfer is not 
completely conservative. We now revisit the possible physical 
foundations for this mass transfer according to the new parame- 
ters of the system derived in previous sections. 

The observations suggest that the mass transfer is driven by a 
wind-Roche Lobe overflow (Podsiadlo wski & Mohamedl2 007). 
We indeed show here that the current state of the system seems 
to require an enhanced mass loss from the giant and that this 
wind possibly fills the Roche lobe and makes the mass transfer 
almost conservative. We also show that it is quite possible that 
an accretion disc formed around the A star, which may explain 
its abnormal luminosity. 



3 IVerh oelst et al] calculate a f actor o f 1.058 . Using the results of 
lHanburv Brown et alj dl974l) and lClaret 



(2000), we find a conversion 



Davis et al ] d2000h . it is equal 



factor of 1.044, and with the m ethod of 1 
to 1.030 ( Ivan Belle et aT]|2009l) . 

4 Since we expect the M star to have had an initial mass at least 
20% higher than the A one q = M M0 /M A0 > 1.2, and therefore 
M M ,o + M A fl > 1.85M M> o. Finally, because the system was most likely 



more massive initially, we get 1.85M MO 
0.55(M M + M A ) =! 2.2 M s . 



> M A + M M , that is, M k 



6. 1 . Mass transfer by Roche Lobe overflow from the wind 

Our results indicate that the M giant only fills around 85 + 
3% of its Roche lobe (Fig. [5). This contradicts the results of 
IVerhoelst et alJ d2007l) . The reasons of this difference are three- 
fold. 

1 . Our more precise interferometric measurements led us to es- 
timate a smaller giant radius than in previous studies. 
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2. The revised HIPPARCOS distance brings the system closer 
than previously thought. This makes the stars smaller, while 
the orbital radius is mostly given by the orbital period and is 
almost independent of the distance. 

3. We have determined for the first time the mass ratio and find 
a higher value than previously guessed. This leads to a higher 
mass and Roche lobe radius of the giant. 

Our results disprove a stricto sensus current RLOF. Two mech- 
anisms might occur that leave this possibility open however, but 
they can be discarded with quantitative arguments. 

First, radiation pressure reduces the gravitational accelera- 
tion influence, so that the Roche pot ential surface shrinks 
( Schuerm anll 19721) . iDermine et al .1 (120091) estimated the ratio of 
the radiation to the gravitational force to be / ~ 10~ 2 to 10 _1 for 
the M giant of SS Lep. It reduces the Roche lobe radius by only 
1 to 4%. For the giant to fill its Roche lobe we need / = 0.35, 
which implies that the L\ and L2 points share the same equipo- 
tentials, making it more difficult for the A star to accrete, most 
of the matter finally going into the circumbina ry disc. However , 
the mass of the latter is rather low according to I Jura et"ail d2001l) 
(Md as t ~ 2 x 1O~ 5 M ), so that this is unlikely. 

Second, atmosphere stratification in red giants indicates that 
there is no single radius v alue, and it may b e not obvious what 
really fills the Roche lobe. lPastetter & Ritterl d!989b have shown 
that for very evolved stars the scale height of the density stratifi- 
cation in their atmosphere is a significant fraction of their photo- 
spheric radius, so that mass flows through the inner Lagrangian 
point L\ long before the photosphere reaches the critical Roche 
equipotential. Because it is on the early AGB phase, the M giant 
in SS Lep is still very far from these evolutionary stages, so that 
this effect is negligible here. 



Finally, to discard the RLOF mechanism, we should ex- 
plain the 0.0126 ma g ellipsoidal variablility in the visible from 
iKoen & Ever! (2002). This variability can be interpreted as the 
signature of a distorted photosphere, whose flux in the line of 
sight varies periodically. In the present case, it would correspond 
to a change in radius of about 8-9 %. However, given the low in- 
clination of the system, even if the star was filling its Roche lobe, 
it would not show su ch a large variation in radius. Moreover, the 
data in the K-band of Kamath & Ashokj (Il999l) . where the M gi- 
ant dominates, are hardly convincing evidence of any periodic 
variations in SS Lep. The Hipparcos variability must thus have 
another origin than a tidal distortion of the M giant. If the reason 
for the light variation in the visible is the primary source of light 
in the system, a change of less than 0.6% is required to explain 
the observed amplitude. The cause of this change is, however, 
not known. 

We therefore conclude that an R LOF is unlikely in SS Lep. 
Podsiadlowski & Mohamed (2007) suggest the possibility of a 
new mode of mass transfer - the wind Roche lobe overflow - 
where a slow wind fills the Roche lobe (e.g. Mira stars in symbi- 
otic systems). Because the wind speed in M giants is rather small 
(around 10 ~ 15kms _1 ) and lower than the orbital one for SS 
Lep (v rb = 48 km s~'), we expect it to be in the particular case of 
a wind Roche lobe overflow, where a substa ntial part of the stel- 
lar wind can be accreted. The simulations of Na gae et alJ (2004) 
show that at least 10% of the M giant wind could be accreted in 
SS Lep. 




Fig. 5. Representation of the modified Roche equipotential (solid 
line) for a mass ratio l/q — 2.2. The limb-darkened diameter of 
the M giant is the dashed line, while the A star one is the dark 
dot. 



6.2. Need for enhanced mass loss from the M giant 

As seen in the previous section, the mass transfer in SS Lep 
is likely due to high-efficiency accretion of the M giant stel- 
lar wind. However, it is difficult to explain the current system 
state with normal stellar wind rates, which are too low. Indeed, 
before the AGB phase, the typical mass loss rates are around 
~l-2 x lO _8 M yr _1 at the normal (e.g. Reimers) rate. The M 
giant (with an expected initial mass > 2.2 M Q ) should have lost 
only a few hundredths of a solar mass before reaching the AGB, 
whereas we expect it to have lost at least 0.9 M . As the M star is 
on the AGB since only a few million years - and will stay there 
for a few million years more at most - it cannot have lost much 
mass since then. 

There is, however, some evidence of enhanced wind mass 
loss of giants i n binaries compared to single gian ts of the same 
spectral type (Mikolaiewska 2007; Jorissen 2003). From a the- 
oretical point of view, the presence of a companion reduces 
the effective gravity of the mass-losing star, thus enhancing the 
mass loss. In the case of SS Lep, the so-called CRAP mecha- 
nism of tidally enhanced stellar wind dTout & Eggletoni ri988) 
allows a mass loss rate 150 tim es higher than the Reimers 
rattfl i.e. ~ 2.4 x 10~ 6 M R yr~^ ISoker. Rappaport. & Harpazl 
( 1998]) and lFrankowski & Tvlendal(l200ll) have also shown that, 
in this case, the mass loss is strongly enhanced in the equatorial 
plane, while an accretion disc can form during win d accretion 
ilTheuns. Boffin.& Jorissenlll996tlNagae et alJl2004 . 

To validate our scenario of a wind RLOF and of en- 
hanced mass loss by wind from the giant, we have consid- 
ered the possible evolution of a binary system, taking the 
CRAP mechanism into account, and following the methodol- 
ogy of iHurlev.Tout. & Polsl (2002) .We start with a system hav- 
ing an initial period of 160 days and initial masses Mm,o = 
2.28 M , M A ,o = 1-85 M Q . For about 1 Gyr, the system evolves 
without much change, and the primary star starts its ascent of 
the AGB. After 2.8 Myr, the masses and period have reached the 
currently observed values, with about 0. 1 M having been lost 
by the system, and forming some circumbinary disc. The mass 
loss and transfer occurred mostly during the last 500 000 years, 
with a mass loss ~ 10~ 6 M yr~'. During the whole process, the 
Roche lobe radius around the M star remained similar, the lowest 
value being 74 R . No RLOF should thus have happened, unless 
the initial eccentricity was very high. 

The above scenario is certainly not the only possible one, but 
it shows that we can explain the current and peculiar properties 
of SS Lep (including a low-mass circumbinary disc), without 



5 Because the M star fills more than 50% of its Roche lobe, we con- 
sider an e nhancement facto r 1 + B/2 6 = 156, where B ~ 10 4 , as de- 
scribed bv lTout & Eggleton] dl988l) . 



6 



N. Blind et al.: An incisive look at the symbiotic star SS Leporis 



resorting to an RLOF and with a low mass loss rate of the order 
of lO^Moyr -1 . 

6.3. Accretion on the A star 

From the previous sections, the accretion efficiency is ex- 
pected to be much higher than 10% in SS Lep, perhaps in the 
rage of 80-90%, while the mass loss rate could reach around 
10 6 M©yr _1 . This corresponds to an accretion rate of a few 
10~ 7 ~ 10~ 6 Mpyr - 1 on the A star This agrees with the 3.3 ~ 
5.5 x 10~ 7 M yr _1 of lBlondel et al] (1 19931) . when interpreting the 
Lyman-a emission in terms of the recombination of H + during 
the inflow of matte r. 

IVerhoelst et all (120071) propose that the A star could have 
"puffed up" by factors of five to seven because of too fast ac- 
cretion in order to explain its abnormal luminosity. The calcu- 
lations of Kippenhahn & Mever-H ofmeisterl (1 19771) require an 
accretion rate of 5 x 10~ 5 M o yr _1 to reconcile the luminosity and 
the spectral type for an increase in diameter of a factor of ten0. 
This is too high compared to the value we derived previously. 
Moreover, in this scenario a shell spectrum is hiding the A star 
spectrum so that a clear determination of the stellar parameters 
is prevented; in particular, the gravity of the star cannot be es- 
timated accurately. If we make the reasonable assumption that 
the A star is rotating i n the orbital plane , based on the measured 
v sin ; = llSkms 1 dRover et al. N 2002). we derive an equato- 
rial velocity of 196kms~ 1 . It is above the break-up speed for a 
2.7 M Q star with a radius of 18 R , which is 170kms _1 . These 
arguments challenge the "puffed up" scenario, so we can wonder 
if the A star indeed has such a large radius. 

Another explanation for the abnormal luminosity of the A 
star could be the presence of an accretion disc. In deed, the Ly-a 
emission profile observed by [Blondel et alj ( 1 1993b consists of a 
single, asymmetric and redshifted feature, suggesting a signif- 
icant absorption by the extended atmosphere in the equatorial 
plane. This could be a hint of wind RLOF presenting a focused 
wind in the equatorial plane and of the presence of an accretion 
disc. This disc would easily explain the shell nature of SS Lep as 
well as the variability of its spectrum. This accretion disc is also 
expected from a mechanical point of view in the wind RLOF 
scenario because most of the matter goes through the L\ point. 
The matter will fall towards the A star and, because of its initial 
angular momentum, it will reach a minimum radius of 

— = 0.0488 q-° M \ 
a 

before forming an accretion disc. In the present case, this 
amounts to 20 R , larger than the star (even if it has expanded). 
Moreover, at the beginning of the mass transfer, r m [ n 9 R Q , 
which is also larger than the radius of an A star on the main se- 
quence (~ 2-3 R ). It is thus likely that a disc has formed. The 
radius of this disc would be about 

— =0.0859 q- 0A2b , 
a 

which is r c i rc 33 R currently for SS Lep. The expected ap- 
parent diameter of this disc is between 0.8 and 1 mas if we take 
the inclination into account. This size is at the limit of detection 
with our VLTI interferometric observations but is within reach 
for CHARA instruments like MIRC or VEGA. We think this al- 
ternative deserves further scrutiny. 

6 Based on the same computations, Ver hoelst et al.1 de- 
rived a value of ~ 2 x 1CT 4 M s yr -1 , but the calculations of 
Kippenhahn & Mever-Hofmeister give an increase in radius by a 
factor of 100 for this accretion rate. 



6.4. Future evolution 

It is also interesting to consider the future of SS Lep in the 
framework of the model presented in Sect. 16. 21 As the M star 
will evolve along the AGB, its mass loss will increase more 
and more. Given the multiplication factor imposed by the CRAP 
mechanism, after 170 000 years, the envelope will be exhausted 
and the star will start its post- AGB evolution to become a WD. 
It is likely that for a short time, the system will appear as a non- 
spherical planetary nebula, where the asymmetry is due to the 
circumbinary material. The A star will have its mass increased 
to 3.3 M . Because the primary will not have time to go through 
the thermal pulses phase, it will most likely not be polluted in 
s-process elements. The period will have increased to 900 days, 
so the system will appear in the typical lo cation in the (e - log P) 
diagr am for post-mass transfer systems dBoffin. Cerf. & Paulusl 
1993). When in the post-AGB phase, the system will thus ap- 
pear typical of those discovered, with a circu mbinary disc and a n 
orbital period between 200 and 1800 days (I van Winckelll2003l) . 
The A star will then undergo its evolution as a red giant and 
the system will most likely again be a symbiotic system - more 
usual this time, given the presence of a white dwarf. 

7. Conclusion and future work 

We have presented here the results of VLTI observations, and we 
focused on the binary. After having computed the characteristics 
of the orbit, we demonstrated that the mass ratio is lower than 
previously thought and that the M giant does not fill its Roche 
lobe. However the system is possibly in the configuration of a 
wind Roche lobe overflow, where a substantial part of the gi- 
ant's stellar wind can be accreted by the A star. We also have 
good reasons to think that the A star is actually surrounded by 
an accretion disc, although this needs to be investigated further. 

We still lack a low-resolution spectro-photometry of SS Lep 
between 1.6 and 2.5 //m to compute the absolute luminosity of 
each component and extract more specific information. The cur- 
rent data also present good hints of unmodelled material escap- 
ing the system. It would be of uttermost importance to obtain 
additional interferometric observations of SS Lep - at large but 
also small baselines - to study the binary interactions in more 
details (M star distortion, mass transfer, circumbinary disc) and 
to determine if the A star is surrounded by an accretion disc. 
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Appendix A: Data and results of the parametric modelling 

The interferometric data are plotted in Figs. lA. lllA.2l and lA.3l (visibilities and closure phases). The reduced data are available on Vizier/CDS. We also plotted on 
these figures the visibilities and closure phases obtained from our binary model with the best parameters. Table I A. H and Fig. lA.4| present the best parameters with 
the data. We indicate for each observation which parameters are free in the fitting procedure. The relative fluxes are always free parameters and are considered to be 
chromatic. Uncertainties on the fitted parameters were computed with Monte-Carlo simulations from the errors on the visibilities and closure phases. 
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Fig. A.l. (u, vj-plane coverage for observations of SS Lep (top: AMBER; bottom: PIONIER). 
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Fig. A.2. Visibility curves for the observations of SS Lep. The order of the figures corresponds to the f«,vj-plan coverages in Fig 
I A. 1 1 In black are the data with error bars. The visibility curves from our model with the best parameters are in red. 
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Fig.A.3. Closure phase curves for the observations for the observations of SS Lep. The order of the figures corresponds to the 
( M,v)-plan coverages in Fig. lA.H In black are the data with error bars. The visibility curves from our model with the best parameters 
are in red. 
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Fig.A.4. Relative flux contribution of the M giant (red), the A star (blue), and the envelope (magenta). The different symbols 
represent the different observations Al (+), A3 (x), PI (□), P2 (A), P3 (V), P4 (O)- 
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Table A.l. Most reliable parameters obtained with the 8 data sets. 



[mas] [mas] [mas] [deg] 

Period Al: m jd ~ 54781 and 54783 (AMBER data) 

Free? Yes Yes Yes Yes 

Value 2.24 12.35 4.2 135 0.71 

Error (1-cr) 2.7 x 1(T 2 2.3 x 1(T 1 3.5 x 1(T 2 0.35 

Period A2: m jd ~ 54826 (AMBER data) 

Free? Yes No Yes Yes 

Value 2.05 12.2 3.8 69 0.20 

Error (l-o-) 2.9 x 10~ 2 - 4.9 x 10~ 2 0.53 

Period A3: m jd ~ 54883 and 54890 (AMBER data) 

Free? Yes Yes Yes Yes 

Value 2.14 12.2 4.5 347 1.36 

Error (l-o-) 4.3 x 10~ 2 2.3 x 10" 1 3.8 x 10" 2 0.37 

Period A4: mjd ~ 54928 (AMBER data) 

Free? Yes No Yes Yes 

Value 1.98 12.2 4.2 291 1.62 

Error (1-cr) 3.9 x 10~ 2 - 4.0 x 10~ 2 0.47 

Period PI: mjd ~ 55500 (PIONIER data) 

Free? Yes No Yes Yes 

Value 2.17 12.2 3.9 32 1.58 

Error (1-cr) 1.4 x 10~ 2 - 1.5 x 10" 2 0.06 

Period P2: mjd ~ 55529 (PIONIER data) 

Free? Yes No Yes Yes 

Value 2.27 12.2 4.5 351 2.46 

Error (l-o-) 1.2 x 10" 2 - 0.2 x 10" 2 0.78 

Period P3: m jd ~ 55537 (PIONIER data) 

Free? Yes No Yes Yes 

Value 2.19 12.2 4.5 346 0.96 

Error (l-o-) 1.6 x 10" 2 - 0.6 x 10" 2 0.20 

Period P4: m jd ~ 55552 (PIONIER data) 

Free? Yes No Yes Yes 

Value 2.30 12.2 4.42 330.6 0.96 

Error (l-o-) 0.9 x 10~ 2 - 0.9 x 10~ 2 0.12 



Notes. The errors are the result of Monte Carlo computations based on the error on the visbilities and closure phases measurements.^ and <f> D 
are the M star diameter and the dusty disc/envelope FWHM, respectively, a and 6 are the visual separation and orientation of the binary. The last 
column gives the final reduced ;f 2 of the fit. 
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